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Abstract 

• We have investigated the spectra of light-heavy quarkonia with the use of a 

quantum-chromodynamic potential model which is similar to that used earlier 
| | for the heavy quarkonia. An essential feature of our treatment is the inclusion 

| of the one-loop radiative corrections to the quark-antiquark potential, which 

m ■ 

t^J- ■ contribute significantly to the spin-splittings among the quarkonium energy 

o\ : 

. levels. Unlike cc and bb, the potential for a light-heavy system has a com- 

"xl" ■ 

plicated dependence on the light and heavy quark masses m and M, and it 

43 ' 

Q_i' contains a spin-orbit mixing term. We have obtained excellent results for the 

observed energy levels of D°, D s , B°, and B s , and we are able to provide pre- 



43 



dieted results for many unobserved energy levels. Our potential parameters 



for different quarkonia satisfy the constraints of quantum chromodynamics. 

We have also used our investigation to test the accuracy of the heavy quark 
effective theory. We find that the heavy quark expansion yields generally good 
results for the B° and B s energy levels provided that M and M hxM 
corrections are taken into account in the quark-antiquark interactions. It does 
not, however, provide equally good results for the energy levels of D° and D s , 
which indicates that the effective theory can be applied more accurately to 
the b quark than the c quark. 
14.40.Lb,14.40.Nd,12.39.Pn,12.39.Hg 
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I. INTRODUCTION 



The light-heavy quarkonia D, D s , B, and B s are of much experimental and theoretical 
interest, and their exploration is necessary for our understanding of the strong as well as the 
electroweak interactions [1,2]. We shall here investigate the spectra of light-heavy quarkonia 
with the use of a quantum chromodynamic model similar to the highly successful model used 
earlier for the heavy quarkonia cc and bb [3]. The complexity of the model is necessarily 
enhanced for a light-heavy system because the potential has a complicated dependence on 
the light and heavy quark masses m and M, and it contains a spin-orbit mixing term. We 
shall obtain results for the observed and the unobserved energy levels of D°(cu), D s (cs), 
B°(bd), and B s (bs), compare them with the available experimental data, and examine their 
scaling behavior. The u-d mass difference and the electromagnetic interaction will be ignored 
in the present investigation [4]. 

We shall also use our results to test the accuracy of the heavy quark effective theory 
[5,6] both in the limit of M — > oo as well as with the inclusion of the M" 1 and M^ 1 InM 
corrections. 

The approximate heavy quark symmetry, like the approximate chiral symmetry, points 
to an underlying difference between the (u, d, s) and the (c, b, t) quarks. It is interesting that 
this fundamental difference was recognized in our mass-matrix approach to quark mixing 
and CP violation [7], which predicted the value M t < 170 GeV for the top quark mass in 
excellent agreement with the recently reported experimental value of 174 ± 17 GeV [8]. 

II. LIGHT-HEAVY QUARKONIUM SPECTRA 

Our treatment for the light-heavy quarkonia is similar to that for cc and bb [3] except 
for the complications arising from the difference in the quark and antiquark masses. Thus, 
our model is based on the Hamiltonian 

H = H + V p + V c , (1) 
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where 

i/ = ( m 2 + p 2 ) l/2 + (M2 + p2) l/ 2 (2) 

is the relativistic kinetic energy term, and V p and V c are nonsingular quasistatic perturba- 
tive and confining potentials, which are given in Appendix A. Our trial wave function for 
obtaining the quarkonium energy levels and wave functions is of the same form as in the 
earlier investigations. Since our potentials are nonsingular, we are able to avoid the use of 
an illegitimate perturbative treatment. 

The experimental and theoretical results for the energy levels of the light-heavy quarko- 
nia D°, B°, D s , and B s , together with the z P[- l P[ mixing angles arising from the spin-orbit 
mixing terms, are given in Tables I-IV. For experimental data we have relied on the Particle 
Data Group [9] except that we have used the more recent results from the CLEO collabo- 
ration [10,11] for D®, Dlf, and D s2 and from the CDF collaboration [12] for B s . In these 
tables, one set of theoretical results corresponds to the direct use of our model, while the 
other two sets are obtained by means of heavy quark expansions of our potentials to test the 
accuracy of the heavy quark effective theory with the inclusion of the M~ x and M~ x InM 
corrections as well as without these corrections. The approximate potentials corresponding 
to the effective theory are given in Appendix B. 

We expect the dynamics of a light-heavy system to be primarily dependent on the light 
quark. Therefore, our potential parameters for D° and B° are the same except for the 
difference in the c and b quark masses, and they are given by 

m ud = 0.350 GeV, 
M c = 1.690 GeV, 
M b = 5.400 GeV, 
fx = 0.932 GeV, (3) 
a s = 0.3965, 
A = 0.185 GeV 2 , 
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B = 0.152. 



Similarly, the parameters for D s and B s are 



m s = 


0.514 GeV, 


M c = 


1.578 GeV, 


Mi, 


5.040 GeV, 


n = 


1.250 GeV, 


a s = 


0.340, 


A = 


0.198 GeV 2 , 


B = 


0.131. 



(4) 



We have ensured that the values of a s , M c , and M b in (3) and (4) are related through 
the quantum chromodynamic transformation relations 

a ' s = l + /5o(a s /4vr)ln(^ 2 //i 2 ) ' (5) 

and 

/ ,x 2 70 /A, 

M' = M i^j , (6) 

where (3 — 11 — |n/, rif = 3, and 7 = 2. The use of the one-loop transformation relations 
is consistent with the inclusion of the one-loop radiative corrections in the quarkonium 
potentials. Moreover, since u, d, and s are the dynamical quarks in the light-heavy systems, 
a higher value of /x for quarkonia with the s quark is to be expected. 

A precise determination of the potential parameters for the light-heavy quarkonia is dif- 
ficult because of the availability of only limited experimental data. This difficulty, however, 
has been mitigated in our treatment by requiring that the parameters for the four systems 
satisfy reasonable physical and quantum-chromodynamic constraints. 

We have also looked at the correlation of our parameters for the light-heavy quarkonia 
with those of other quarkonia. When applied to ud and us, our parameters yield good results 



for the ir-p and K-K* splittings. We are, however, unable to correlate our parameters for 
the light-heavy quarkonia with those for the heavy quarkonia through the transformation 
relations (5) and (6), and keeping in mind the past success of our quarkonium model we can 
only offer the following possible explanation: 

Strictly speaking, the QCD transformation relations are applicable only to the current 
quarks. According to our experience, the transformation relations seem to hold reasonably 
well for the heavy quarkonia cc and bb as well as for the quarkonia containing one or two 
light quarks, but they are unsuitable for correlating the parameters for these two classes 
of quarkonia. This seems to be another manifestation of the difference between light and 
heavy quarks. We believe a full explanation would require an understanding of the origin of 
the constituent quark masses, which remains unclear at this time. 

Our phenomenological confining potential for the light-heavy quarkonia is of the same 
form as that for the heavy quarkonia. We find that the parameter A for the spin-independent 
term in the confining potential is approximately the same for all quarkonia, while the spin- 
dependent terms vary such that the vector-exchange component is smaller for cc than for 
bb, and still smaller for the light-heavy quarkonia. 

III. CONCLUSION 

We have obtained excellent results for the observed energy levels of D°, B°, D s , and 
B s with the use of our quantum-chromodynamic potential model, and provided predicted 
results for many unobserved energy levels in Tables I- IV. We have included in these tables 
the mixing angles for the 1 3 P[ and 1 1 P[ levels, which are needed for an understanding of 
their decay properties. Although the use of a semirelativistic model may seem questionable 
for a system containing a light quark, ultimately such an approach should be judged on 
the basis of its predictions [13]. Additional experimental data on the light-heavy quarkonia 
should be available in the near future. 

We have also used our results to test the accuracy of the heavy quark effective theory [5,6]. 
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By comparing the theoretical results without and with the effective theory in Tables I- IV, we 
find that the heavy quark expansion with the inclusion of the M _1 and M _1 In M corrections 
yields generally good results for the B° and B s energy levels. It does not, however, provide 
equally good results for the energy levels of D° and D s , which indicates that the effective 
theory can be applied more accurately to the b quark than the c quark [14]. 

We further find that the results for the energy levels in the limit M — > oo are unac- 
ceptable. As is well known, in this limit the energy level pairs ( 3 S'i, 1 S'o), ( 3 P[, 3 Po), and 
( 3 P 2) 1 P[) become degenerate. 

Finally, we have examined the scaling behavior of energy level splittings in the light- 
heavy quarkonia by looking at the results obtained by the direct use of our model in Tables 
I-IV. As shown in Tables V and VI, the splittings between levels which become degenerate 
in the limit M — > oo exhibit an approximate M _1 scaling. This scaling behavior does not 
apply to splittings between other pairs of energy levels. 
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APPENDIX A: NONSINGULAR QUARKONIUM POTENTIALS 

The nonsingular potentials for a light-heavy quarkonium are similar to those used recently 
for cc [3] except for the complications due to the difference in the quark and antiquark masses. 
The complications are further enhanced by the conversion of the singular potentials [15] into 
the nonsingular ones [16], which are necessary to avoid the use of an illegitimate perturbative 
treatment. The corresponding denominators in the singular and nonsingular potentials for 
a quark and an antiquark of masses mi and m 2 are related as 

11 11 



mf mf + |k 2 ' m| m| + 7k 2 
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mi m 2 ~* {ml + ik 2 )V2( m 2 + l k 2)l/2 ~ mim2 + 1 frrn + k 2 

o y 7712 7Ti-]_ y 



or 



(A2) 



171x1712 + |k 2 

where 

k = — o ~ , a; = {K,mim 2 ) 1 . (A3) 

mf + mi 

The potentials for a quark and an antiquark of different flavors are given below. They 
reduce to those for cc for m 1 = m 2 = m except that, unlike cc, they do not contain the 
annihilation terms. 
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In the coordinate space, it takes the form 
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Note that the tensor operator is defined as 
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the functions E± are expressible in terms of the exponential-integral function Ei as 
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2. Phenomenological confining potential 



The scalar-vector-exchange confining potential is given by 
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The coordinate-space potentials are 
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It is understood that the confining potential also contains an additive phenomenological 
constant C . 



APPENDIX B: QUARKONIUM POTENTIALS WITH HEAVY QUARK 

EXPANSION 

Upon replacing mi and m 2 by m and M, and expanding in powers of M~ l , the coordinate- 
space potentials of Appendix A take the approximate forms given below. 

1. Perturbative Potential 
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2. Confining potentials 
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3. M -»■ oo limit 

In the limit M — > oo, the perturbative and confining potentials are given by 

V p = V p0 , V S = V S0 , V V = V V0 . (B9) 

According to (B2), (B5), and (B7), the spin-orbit and spin-orbit-mixing terms in these 
potentials are of the form 

/(r)L.S + /(r)L.(S!-S 2 ), (BIO) 

which is expressible solely in terms of the light-quark spin as 

2/(r)L-S 1 . (Bll) 

We have, however, given the potentials in the form (BIO) to facilitate comparison with the 
more accurate treatments of the light-heavy quarkonia. 
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TABLES 

TABLE I. D° energy levels in MeV. Effective theory results are given with the M _1 and 
M _1 InM corrections as well as in the limit of M — ► oo. Experimental results are from Refs. [9] 
and [10]. 
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TABLE II. D s energy levels in MeV. Experimental results are from Refs. [9] and [11]. 





Expt. 


Theory 


Effective theory 


M -> oo 


i % (AO 


1968.8±0.7 


1968.8 


1968.8 


1968.8 




2110.3±2.0 


2110.5 


2113.1 


1968.8 


2 % 




2656.5 


2678.8 


2536.5 


2 3 Si 




2757.8 


2774.3 


2536.5 


1 13 P 




2387.8 


24:22.2 


o o o o o 

2382.2 


1 3 P{ 




2521.2 


2528.8 


2382.2 


1 3 P 2 (£> s2 ) 


2573.2±1.9 


2573.1 


2582.8 


2402.8 


1 ^ (Al) 


2536.5±0.8 


2536.5 


2552.1 


2402.8 


9 




26.0° 


31.8° 


35.6° 


TABLE III. 


B° energy levels 


in MeV. Experimental results are from Ref. [9]. 






Expt. 


Theory 


Effective theory 


M -> oo 


1% (5°) 


5278.7±2.1 


5278.7 


5278.7 


5278.7 


1 3 6"i (5*°) 


5324.6±2.1 


5324.0 


5325.8 


5278.7 


2 




5892.1 


5893.9 


5846.3 


2 3 Si 




5924.3 


5927.1 


5846.3 


1 3 P 




5689.5 


5692.5 


5659.1 


1 3 P{ 




5730.8 


5734.1 


5659.1 


1 3 P 2 




5759.1 


5761.4 


5701.5 


1 !P{ 




5743.6 


5745.4 


5701.5 







31.7° 


31.3° 


35.6° 
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TABLE IV. B s energy levels in MeV. Experimental results are from Refs. [9] and [12]. 





Expt. 


Theory 


Effective theory 


M -> oo 


1 % (B s ) 


5383.3±6.7 


5383.3 


5383.3 


5383.3 


l 3 Si {Bt) 


5430. 5±2. 6 


5431.9 


5434.1 


5383.3 


2 1 S 




6000.9 


6003.1 


5950.9 


2 3 Si 




6035.8 


6039.1 


5950.9 


1 J P 




5810.1 


5814.2 


5796.7 


1 3 P{ 




5855.0 


5857.9 


5796.7 


1 3 P 2 




5875.2 


5878.1 


5817.1 


1 !p{ 




5860.2 


5863.2 


5817.1 






27.3° 


27.1° 


35.6° 


TABLE V. 


Scaling of energy level 


splittings in D° 


and B°. Splittings are given 


in MeV. 




AL>o 




(m-c/mj ) lAJJ 




1 3 5i - 1 % 


142.5 




44.6 


45.3 


2 3 5i - 2 % 


99.3 




31.1 


32.2 


1 3 P[ - 1 3 P 


128.7 




40.3 


41.3 


1 3 P 2 - 1 l P[ 


44.0 




13.8 


15.5 


TABLE VI. Scaling of energy level 


splittings in D s 


and B s . Splittings are given 


in MeV. 




AD S 




(m c /m b )AD s 


A T~> 

AB S 


1 3 5i - 1 l S 


141.7 




44.4 


48.6 


2 3 Si - 2 % 


101.3 




31.7 


35.0 


1 3 P{ - 1 3 P 


133.4 




41.8 


44.9 


1 3 P 2 - 1 *P{ 


36.6 




11.5 


15.0 
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